arXiv:astro-ph/0503063vl 3 Mar 2005 


To appear in ApJ 


Multi-wavelength carbon recombination line observations with 
the VLA toward an UCHll region in W48: Physical properties 
and kinematics of neutral material 

D. Anish Roshi 

Raman Research Institute, Sadashivanagar, Bangalore, India 560080, anish@rri.res.in 

and 

National Radio Astronomy Observatory 1 , Green Bank, USA 

W. M. Goss 

National Radio Astronomy Observatory, P.O. Box 0, Socorro, NM 87801, USA; 

mgoss@nrao.edu 

K. R. Anantharamaiah 2 

Raman Research Institute, Sadashivanagar, Bangalore, India 560080 

and 

S. Jeyaknmar 3 

Raman Research Institute, Sadashivanagar, Bangalore, India 560080 

ABSTRACT 

Using the Very Large Array (VLA) the C76a and C53a recombination lines 
(RLs) have been detected toward the nltra-compact Hu region (UCHll region) 
G35.20—1.74. We also obtained upper limits to the carbon RLs at 6 cm (CllOa & 
Cilia) and 3.6 cm (C92a) wavelengths with the VLA. In addition, continuum 
images of the W48A complex ( which includes G35.20—1.74 ) are made with 
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angular resolutions in the range 14" to 2". Modeling the multi-wavelength line 
and continuum data has provided the physical properties of the UCHll region 
and the photodissociation region (PDR) responsible for the carbon RL emission. 
The gas pressure in the PDR, estimated using the derived physical properties, 
is at least four times larger than that in the UCHll region. The dominance of 
stimulated emission of carbon RLs near 2 cm, as implied by our models, is used 
to study the relative motion of the PDR with respect to the molecular cloud 
and ionized gas. Our results from the kinematical study are consistent with a 
pressure-confined UCH II region with the ionizing star moving with respect to the 
molecular cloud. However, based on the existing data, other models to explain 
the extended lifetime and morphology of UCHll regions cannot be ruled out. 

Subject headings: ISM: Hll regions — ISM: general - radio continuum: ISM - 
radio lines: ISM - line: formation 


1. Introduction 

Massive stars (OB stars) are formed by the gravitational collapse of clumps in molecular 
clouds. Lyman continuum (E > 13.6 eV) photons from newly born OB stars ionize the sur¬ 
rounding material resulting in the formation of ultra-compact Hll regions (UCHll regions). 
LICH II regions are expected to expand due to their high gas pressure. If they expand at the 
speed of sound in the ionized gas, then the time taken for the UCHll region to expand to 
~ 0.1 pc is a few times 10 3 years. This time scale is referred to as the dynamical lifetime. 
However, the lifetime deduced from the observed number of UCH II regions is a few times 10 5 
years. This discrepancy between the two lifetimes is referred to as the “lifetime problem” of 
UCHll regions (Wood & Churchwell 1989a). Resolving the lifetime problem is an important 
step in understanding massive star formation in the Galaxy (see review by Garay & Lizano 
1999). 

There have been several suggestions to resolve the lifetime problem. In particular, De 
Pree, Rodriguez & Goss (1995) proposed that if high density (~ 10' cm -3 ), warm (~ 100 
K) molecular material is present in the vicinity of UCHll regions, it may be able to pres¬ 
sure confine LTCHlI regions that form there and thus extend their lifetime. In the recent 
past, attempts have been made to observationally determine the physical properties of the 
molecular material near UCHll regions. For example, Akeson & Carlstrom (1996) have used 
methyl cyanide, a good tracer of the temperature and density of dense molecular cores, to 
estimate the physical properties of the ambient medium near UCHll regions G5.89+0.4 and 
G34.3+0.2. They concluded that the estimated ambient pressure (> 10 8 K cm -3 ) was high 
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enough to pressure confine the UCHlI regions. The dominant component of the ambient 
pressure need not always be the thermal pressure of the molecular gas but can be the tur¬ 
bulent pressure (Xie et ah 1996). Advances have also been made in modeling the dynamical 
evolution of UCH II regions in the presence of high pressure molecular environment. Through 
simulations, Garcia-Segura & Franco (2004) show that pressure-confined UCH II regions with 
ionizing stars moving with respect to the molecular core can be long-lived and can produce 
the observed morphologies of the H II regions. 

If UCH II regions are pressure confined due to high density gas in their vicinity, then 
other observable effects will occur. In particular, far-ultraviolet (FUV) photons (6.0 - 13.6 
eV) from the OB star would produce photo-dissociation regions (PDRs) in the neutral ma¬ 
terial close to the UCH II region. Near the interface between the region where hydrogen is 
ionized and the PDR, gas phase carbon will be ionized by FUV photons in the energy range 
11.3 to 13.6 eV (see review by Hollenbach & Ticlens 1997). The electron density in this 
layer is relatively high (> 10 3 cm -3 ) and gas temperatures can be in the range 300 - 1000 
K (e.g. Natta, Walmsley & Ticlens 1994). These conditions are ideally suited for producing 
observable radio recombination lines (RLs) of carbon. 

Physical properties of dense material near UCHlI regions were earlier estimated from 
observations of high density molecular tracers. However, high line optical depth in these 
regions complicates the determination of the physical properties. We propose that multi¬ 
wavelength carbon RL observations from PDRs associated with UCHlI regions can be used 
to estimate the properties of the dense molecular material. Unlike molecular traces, RLs do 
not suffer from the limitation of high opacity and thus provide another probe to test whether 
LICHll regions are pressure confined. 

In this paper, we present multi-wavelength (0.7, 2, 3.6 & 6 cm) VLA observations 
of carbon RLs from the PDR associated with G35.20—1.74. The UCHlI region is in the 
molecular cloud complex W48 at l = 35°.2, b = —1°.74 at a distance of 3.2 kpc (Wood & 
Churchwell 1989b). At 0.7 cm, we have detected the C53a and X53a (see §4) transitions. 
This data forms the first successful imaging of carbon RL at 0.7 cm with the VLA at an 
angular resolution of about 2". The details of the observations are given in §2. In §3 and §4, 
we discuss the modeling of the multi-wavelength data to determine the physical properties 
of the PDR and the UCHlI region. In §5, we estimate the gas pressure in the UCHlI region 
and the PDR and investigate whether the HII region is pressure confined. 
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2. Observation and data reduction 

We made spectroscopic observations of W48 using the VLA in D-configuration at 0.7, 2, 
3.6 and 6 cm in dual polarization mode. The 0.7 cm observations were made during August 
2004 and data at other wavelengths were obtained during October 2001. The transitions 
observed are the CllOa (4876.5886 MHz) and Cilia (4746.5497 MHz) at 6 cm, the C92a 
(8313.5279 MHz) at 3.6 cm, the C76a (14697.3141 MHz) at 2 cm, and the C53a (42973.3984 
MHz) at 0.7 cm. The bandwidth at 2 cm wavelength was 3.13 MHz, which corresponds to a 
total velocity coverage of 64 km s _1 . Thus in addition to the carbon line, roughly 75 % of the 
He76a line profile falls within the observed bandwidth. Table 1 summarizes the parameters 
of the observations. 

Data analysis was carried out using the Astronomical Image Processing Software (AIPS). 
The default channel zero data was used for continuum calibration. After satisfactory editing 
and calibration, the flag and calibration tables were transferred to the spectral data. The 
system band shapes were obtained using the bandpass calibrator data with the AIPS task 
BPASS. The estimated band shapes were used for bandpass calibration. Line free channels 
from the bandpass calibrated data were used for estimating the continuum emission which 
was then subtracted from the spectral line data. We used the task UVLSF for this purpose. 
The continuum uv data created by UVLSF was used as the input to IMAGR to make the 
continuum images. The spectral cubes were also made with the AIPS task IMAGR. The 
data at 0.7 cm were analyzed by taking into account of the weighting provided by the on-line 
system for each visibility measurement. GIPSY (Groningen Image Processing System) and 
AIPS++ software packages were used to further process (e.g. Gaussian fit) the spectral line 
data. 


3. Continuum Emission 

Fig. 1 shows continuum images at 0.7, 2, 3.6 and 6 cm with angular resolutions ranging 
from ~ 2" to 14". These continuum images show a compact source of angular size ~ 5" 
x 5" along with an extended source (angular size ~ 1'.5 x F.5), located south and east of 
the compact component (see Fig. 1). This extended source is designated as W48A in the 
literature (Onello et al. 1994). It has been suggested that the extended emission is directly 
associated with the compact component (Kurtz & Franco 2002). 

The continuum emission and the partially ionized gas associated with W48A have been 
extensively studied earlier with a lower angular resolution of ~ 60" (see Onello et al. 1994). 
In our images, the emission from W48A extends to about 2' at the longest wavelength. The 
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flux densities of the extended emission measured are 5, 10 and 12 Jy respectively at 2, 3.6 
and 6 cm wavelengths. The measured flux density at 2 cm is about a factor of two less 
than that estimated using the flux densities at 3.6 and 6 cm wavelengths and considering a 
spectral index of —0.1 for the thermal emission. This lower flux density is due to the missing 
short spacings in interferometric observations. The extended emission is not detected at 0.7 
cm because of the same reason. The largest angular size that our observations are sensitive 
to are 5', 3', 1.5' and 0.7' respectively at 6, 3.6, 2 and 0.7 cm. Therefore at 2 and 0.7 cm 
wavelengths the full extent and flux density of the emission from W48A are not determined. 

The compact source in Fig. 1 is an UCHll region referred to as G35.20—1.74 in the 
literature (e.g. Wood and Churchwell 1989b). The continuum emission from G35.20—1.74 
shows some resemblance to a cometary morphology as inferred earlier by Wood & Churchwell 
(1989b). The flux densities of the compact source at the four observed wavelengths are given 
in Table. 2. The angular size of the source is estimated from the 2 and 0.7 cm images, which 
have the highest angular resolutions. The AIPS task JMF1T was used to determine the 
angular size. The flux density and angular size are consistent with those estimated earlier 
(Wood & Churchwell 1989b, Woodward, Heifer & Pipher 1985). 

The continuum data toward G35.20—1.74 are used to determine the physical properties 
of the UGH II region. We model the continuum emission by considering that the emission 
from the compact source originates from a spherical, homogeneous ionized region. The flux 
density due to thermal emission from such an ionized region is estimated as described by 
Mezger & Henderson (1967). The angular size of the spherical region estimated from the 
observed size of the UCHll region is 6". 3. Table 3 gives the parameters of the ionized gas 
obtained from continuum modeling. A plot of the flux density from the model as a function 
of frequency is shown in Fig. 2. At a distance of 3.2 kpc (Wood & Churchwell 1989b), the 
diameter of the spherical region of ionized gas is 0.1 pc. The derived physical properties 
and the observed flux densities are then used to estimate the excitation parameter, rate of 
Lyman continuum photons required to maintain ionization equilibrium and the spectral type 
of the embedded star (Panagia 1973). These values are also included in Table 3. 


4. Recombination Line Emission toward G35.20—1.74 
4.1. Carbon Recombination Lines 

We detected the C53o; and C76a recombination lines in the direction of G35.20—1.74. 
No carbon RLs were detected toward the extended H II region W48A. The non-detections 
toward other directions indicate that the carbon line originates from the PDR associated 
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with the UCHll region. The spectra at 2 and 0.7 cm wavelengths, averaged over a 6".3 
x 6".3 region near the continuum peak toward G35.20—1.74, are shown in Fig. 3 and the 
parameters obtained from the Gaussian fit to the line features in the figure are given in 
Table. 4. The upper limits on the carbon line emission at 3.6 and 6 cm are included in 
Table. 4. The data at 3.6 cm are affected by a systematic baseline ripple that was not 
removed by a careful bandpass calibration. This ripple is present in the spectra toward 
regions with bright continuum emission. The upper limit we obtained toward the UCHII 
region is from the spectrum with the baseline ripple and therefore its value is higher (5.2 
mJy/beam) compared to that estimated from an off source position (= 1.1 mJy/beam). 


4.2. Models for the Carbon line emission toward G35.20—1.74 

We consider homogeneous ‘slabs’ of PDR material placed in front and back of the 
UGH II region and solved the radiative transfer equation for non-LTE cases to obtain the 
recombination line flux density. The non-LTE departure coefficients b n and f3 n are calculated 
using the program originally developed by Brocklchurst & Salem (1977) and later modified 
by Walmsley & Watson (1982) and Payne, Anantharamaiah & Erickson (1994). The two 
relevant input parameters for the program are: (1) abundance of gas phase carbon and 
(2) background radiation held. The abundance of carbon is taken as 0.75 of the standard 
abundance (3.9 x 10” 4 ; Morton 1974), which implies a depletion factor of 25 % (Natta et 
al. 1994). The departure coefficients depend on the background radiation held. We have 
selected a thermal background due to an UGH II region with temperature and emission 
measure determined from the continuum observations (see §3). The departure coefficients 
are computed for a set of electron temperatures in the range 100 to 1000 K and densities 
between 1000 and 6000 cm -3 . The dielectronic like recombination process that modifies 
the level population (Walmsley & Watson 1982) is also included in the calculation. The 
coefficients are computed by considering a 10000 level atom with the boundary condition 
b n —► 1 at higher quantum states (see Payne et al. 1994 for further details). 

The line intensity is a function of PDR gas temperature, Tpdr, electron density, n RDR , 
PDR thickness along the line of sight, l, and the background radiation held. For the homo¬ 
geneous PDR, the line brightness temperature, Tlb, due to the slab in the near side of the 
UGH II region is given by (Shaver 1975) 

Tlb = Tbg, u + T inv , (1) 

where 

T bg , u = T ob g tU e~ TCl '(e~ TLl/ — 1 ) , 
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T in , u = T pdr f hmT ^ + T ^ (l _ e-fa.+rc,)) _ (1 _ e -rc^\ (2) 

V T Lu + T C U J 

Ty XJ)V is the contribution to the line temperature clue to the background radiation held and 
T m i , is the intrinsic line emission from the slab. In Eq. 2, T {)bg u is the background radiation 
temperature, which, as discussed above, is the continuum emission from the UCHll region. 
Tcv is the continuum optical depth of the PDR. The non-LTE line optical depth of the 
spectral transition from energy state m to n, Tl„, is tl v = b n /3 n rl u where t Pv is the LTE 
line optical depth, b n and are the departure coefficients of state n. For the PDR, t Pv oc 
n RDR nc+l, where ric+ is the carbon ion volume number density in the PDR. For the present 
calculations we assumed that n RDR = nc+ = n e so t Rv oc n 2 e l. The line temperature from 
the PDR on the far side (see below) is obtained from Eq. 1 by setting T b(JV = 0. The line 
brightness temperature is finally converted to flux density using an angular size of about 
3". 7 (see Table 2). 

We considered three classes of models - (a) models with line emission due to a PDR slab 
in the front of the UCHll region; (b) models with line emission having contribution from 
PDR slabs in the front and back of the UCHll region and (c) model with line emission due 
to a PDR slab in the far side of the UCHll region. The line emission for class (a) models 
are obtained using Eqs. 1 and 2. For class (b) models, the intrinsic contribution (ie T in ^ 
in Eq. 1) from the PDR slab on the far side is added to Eq. 1. Note that in this case we 
essentially assume that the relative motion between the PDR slabs in the front and back 
of the UCHll region is smaller than the width of the observed carbon line (see §5.2). For 
each class of models, a set of parameter values is determined such that the observed line 
flux densities at 0.7 and 2 cm (also consistent with the observed upper limits at 3.6 and 6 
cm) are reproduced within ± la error. We found that class (c) models cannot reproduce the 
observed line flux densities and hence will not be discussed further. Table 5 gives a subset 
of model parameters that are consistent with our RL data at the four wavelengths. The 
neutral densities listed in Table 5 are obtained with a gas phase carbon abundance of 3 x 
10 -4 used for modeling. Fig. 4 shows the model carbon line flux density as a function of 
frequency along with the observed values and limits for class (a) and (b) models with typical 
temperatures 200 and 500 K and electron density 2500 cm -3 . 

The modeling shows that: (1) the electron density and PDR thickness are well con¬ 
strained by our RL data with the density in the range 1500 - 6000 cm -3 and the PDR 
thickness a few times 10~ 4 pc; (2) models with temperatures < 150 K are ruled out by our 
RL data. We have explored models with PDR temperatures up to 1000 K (Natta et al. 1994). 
However, the temperature could not be well constrained due to the uncertainty in the upper 
limit obtained from the 3.6 cm spectrum (see §4.1). Note that models with temperature < 
250 K are consistent with the 3cr upper limit obtained from the off source region of 3.3 mJy 



at 3.6 cm (see Fig. 4). 


In a recent single dish survey of carbon lines near 8.5 GHz, Roshi et al. (2005) de¬ 
tected RLs toward a large number of UCHll regions. These lines along with upper limits at 
other frequencies were used to model the properties of the line forming region. The ranges 
of density and size obtained for the PDR toward G35.2—1.74 are comparable with those 
determined by Roshi et al. (2005) toward other UCHll regions. 

An important inference from modeling is that the line emission at frequencies near ~ 14 
GHz is dominated by stimulated emission due to the background continuum emission arising 
from the UCHll region. The intrinsic line flux density from the slab near 14 GHz is < 20 
% of the observed RL flux density. Thus for these frequencies Eq. 1 can be approximated as 
Tlb rsj Tbg tU . At frequencies >14 GHz, the intrinsic line emission from the slab (ie T in u in 
Eq. 1) contributes significantly. For example, at 0.7 cm the intrinsic line flux density from 
the PDR slab is about 6 mJy for models with electron density 3000 enr 3 , which is almost 
equal to the line flux density due to the background term (ie Tf )giJ ) in Eq. 1. 


4.3. Helium and Other Recombination Lines 

In addition to the C76a transition, He76a RLs in the direction of G35.20—1.74 and 
W48A were detected. The helium line parameters are given in Table 4. The LSR velocity 
of the He76a line (+42.0±0.5 km s _1 ) differs by —5.9 km s -1 from the earlier measured 
value for the H76o line (Wood & Churchwell 1989b). The parameters obtained for the H76a 
RL are Vlsr = +47.9+1.2 km s -1 , AV = 31.8±1.89 km s^ 1 and Sl = 340±26 mJy. The 
angular resolution of our 2 cm observation is similar to that of the earlier observations (~ 
4"; Wood & Churchwell 1989b). Differences in the sizes of the regions where helium and 
hydrogen are ionized could be the cause of the LSR velocity difference between the two RLs. 
The limited velocity coverage of the 2 cm data (see §2) could also be the cause of the LSR 
velocity difference. Further confirmation is needed however. 

The line emission from W48A is averaged over an area of 36" (in RA) x 17" (in Dec) 
centered at RA(2000) 19 h 01 m 46 s .4 and DEC(2000) 01°13 , 02”. This region does not include 
the UCHll region and thus is an estimate of the helium line emission from W48A. The line 
parameters given in Table 4 are from the average spectrum. The line flux density from W48A 
is only about 6 % of that toward the UGH II region. The LSR velocity of the helium line from 
W48A (+47.4+0.5 km s -1 ) is similar to that of the H76a line detected from G35.20—1.74 
by Wood & Churchwell (1989b). 

At 0.7 cm wavelength, a line feature at LSR velocity about +43.4 km s' 1 is detected. 
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The LSR velocity is obtained by assuming that the mass of the ion to be infinity. Valee (1987) 
have detected the S125a recombination line toward W48A. However, the LSR velocity of 
the sulfur line (+43.1 km s” 1 ; Valee 1987) is similar to that of carbon lines detected in our 
observations. Therefore the second line feature detected at 0.7 cm wavelength may not be 
a sulfur line. The possibility of this line being a Doppler shifted carbon RL cannot be ruled 
out. Further investigation is needed to identify this line feature. 


5. Is G35.20—1.74 pressure confined ? 

As discussed in §1 the lifetime of UCHll regions can be extended to a few times 10 5 
years if they are pressure confined. In this section, the derived physical properties of gas 
inside and outside G35.20—1.74 are used to investigate whether the UCHll region is in fact 
pressure confined. 


5.1. Estimation of gas pressure 

The total gas pressure inside the UCHll region is the sum of thermal pressure and 
turbulent pressure. No measurement of the magnetic field in G35.20—1.74 or the PDR 
associated with the Hll region exists and hence we do not include its contribution in the 
calculation of total pressure. The total gas pressure is given by 

Pucka = 2kn e T e + n e nm H v 2 H _ turb dyne cm” 2 , (3) 

where k is the Boltzmann constant, n e is the electron density in cm -3 , T e is the electron 
temperature in K, ra# is the hydrogen mass in gm and Vu-turb is the turbulent velocity 
inside the LICHlI region in units of cm s” 1 . The line profile due to turbulence is considered 
to be Gaussian with a = vn-turb, which is estimated from the observed FWHM of H76a 
transition (Wood & Churchwell 1989b) and the calculated line width due to thermal motion. 
The effective mass in amu of H + He gas with He fraction in number of atoms taken as 
10 % of that of H is 1.4. From the flux densities of H76a and He76a RLs we infer that 
about 8 % of the helium is ionized. Thus the parameter /r in the turbulent pressure term 
in Eq. 3 is ~ 1.3 since it is expressed in terms of n e . The values for electron density and 
temperature, estimated from modeling the continuum emission from the LICHlI region, are 
used to estimate gas pressure. The total pressure inside the UGH II region is 1.3 x 10”' dyne 
cm” 2 , with a contribution of 7.0 x 10” 8 dyne cm” 2 from thermal process and 5.7 x 10” 8 
dyne cm” 2 from turbulence. The turbulent pressure is about 85 % of the thermal pressure. 
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The total gas pressure in the PDR is given by 

Ppdr = kn H T PDR + n H iim H v 2 c _ turb dyne cm ' 2 , (4) 

where rip is the number density of hydrogen atoms in cm ' 3 , Tpp> R is the PDR gas temper¬ 
ature in K, Vc—turb is the turbulent velocity in cm s' 1 . The line profile due to turbulence is 
considered to be Gaussian with a = vc-turb, which is estimated from the observed FWHM 
of C76 a transition. The effective mass p in amu is taken as 1.4. To estimate PDR pressure 
we consider that in the region where carbon is ionized, the neutral material is predominantly 
in atomic hydrogen form and the gas temperature is equal to the electron temperature es¬ 
timated from carbon line modeling (see Hollcnbach & Teilens 1997). The estimated total 
gas pressure in the PDR is between 5.3 x 10“' and 4.3 x 10 “ 6 dyne cm ' 2 for the model 
parameters given in Table 5. For PDR temperatures < 500 K, the turbulent pressure is at 
least 10% more than the thermal pressure. This importance of turbulent pressure was noted 
earlier by Xie et al. (1996). 

Comparing gas pressures in the UCHlI region and the PDR indicates that the pressure 
in the PDR is at least four times larger than that in the UGH II region. In §5.2, we investigate 
the relative motion between the different components (ionized gas, material in the PDR and 
molecular gas) along the line-of-sight to understand this pressure difference. 


5.2. Kinematics 

The UGH II region G35.20—1.74 is located in the molecular cloud complex W48. The 
carbon line emission originates from PDR formed at the interface between the molecular 
cloud and the H II region. The PDR resides within the molecular cloud. The relative motion 
along the line-of-sight between the three components (ionized gas, material in the PDR 
and molecular gas) can be determined by examining the central velocities of tracers of the 
different components (see Fig. 5). The central velocities of high density tracers of molecular 
gas observed toward G35.20—1.74 and the angular resolution of these observations are given 
in Table 6 . The measured velocities with respect to the Local Standard of Rest (LSR) 
ranges from +42.8 to +46.1 km s' 1 . (The errors in the central velocities are not given in the 
literature except for NH 3 ( 2 , 2 ), which is 0.1 km s' 1 .) A variety of causes such as turbulence, 
shocks, complex excitation and chemistry in the vicinity of the UGH II region contribute to 
the spread in the central velocity. Moreover, the different angular resolutions of the molecular 
line observations may also contribute to the velocity difference. Molecular line observations 
with higher angular resolution may help to precisely determine the LSR velocity of the cloud 
associated with the UCHlI region. Here we take the mean value (+44.4+1.0 km s' 1 ) of the 
central velocities listed in Table 6 as the LSR velocity of the molecular cloud. 
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The H76 a line (Wood & Churchwell 1989b), which originates from the ionized gas in the 
UCHll region (see §4.3), has a LSR velocity of +47.9±1.2 km s _1 . Since the UCHll region 
is optically thin at 2 cm, this velocity provides the mean velocity of the ionized gas with 
respect to the LSR. As discussed in §4.2, the carbon RL flux density at 2 cm is dominated 
by stimulated emission and hence we preferentially observe the emission from the PDR in 
front (ie toward the observer; see Fig. 5) of the LICHll region. The carbon RL at 2 cm has 
a central velocity of +41.9±0.4 km s _1 . 

Translating the measured radial velocities to a reference frame at rest with respect to the 
molecular cloud gives the picture of the PDR moving into the molecular cloud at — 2.5±1.1 
km s -1 and the ionized gas moving at +3.5±1.6 km s^ 1 relative to the molecular cloud (see 
Fig. 5). 


5.3. Discussion 

Modeling the multi-wavelength (A = 0.7, 2, 3.6 & 6 cm) carbon recombination line and 
continuum data toward the UCHll region G35.20—1.74 has provided the physical properties 
of the ionized gas and the PDR associated with the H II region. These physical properties 
were then used to determine the pressure in the LICHll region and the PDR in order to 
investigate whether the UCHll region is pressure confined. Our calculations have shown 
that the pressure in the PDR is at least four times larger than that in the UCHll region. 
Does this large pressure in the PDR indicate that the UCHll region is pressure confined 
? As shown in §5.2, the PDR moves relative to the molecular cloud at —2.5±1.1 km s -1 . 
This relative motion is inconsistent with a ‘stationary’ pressure-confined nebula where the 
H II region is confined within the molecular cloud and the massive star ionizing the nebula 
is stationary with respect to the molecular cloud. 

Garcia-Segura & Franco (2004) have considered the effects due to the natal molecular 
cloud gravity and stellar motion in their recent gas dynamical simulations. These simula¬ 
tions were done to study the evolution of pressure-confined H II regions. They show that 
a plethora of structures for the ionized gas can be produced due to the stellar motion. In 
such models relative motions between the PDR, ionized gas and molecular cloud are ex¬ 
pected both clue to stellar motion and re-adjustments necessary for the ionized gas to be in 
hydrostatic equilibrium in the presence of cloud gravity. The PDR may reside in the shock 
formed as a result of stellar motion, which also explains the high pressure calculated in the 
PDR compared to the ionized gas. If the extended Hll region W48A is directly associated 
with G35.20—1.74 (Kurtz & Franco 2002), then, in this picture, the star has already moved 
into a density ramp in the molecular cloud which has resulted in the formation of the blister 
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type region (Garcia-Segura & Franco 2004). However, other models that are proposed to 
explain the extended lifetime of UCHll regions cannot be ruled out using the present data. 
For example, in the model by Kim & Koo (2001), where they combine the champagne flow 
model with the hierarchical structure of molecular cloud, relative motion of ionized gas with 
respect to the molecular cloud is expected due to the flow. It is possible that the PDR 
is moving with respect to the molecular cloud either because of stellar motion or due to 
increased pressure inside the H II region caused by a stellar wind. 


We are grateful to the anonymous referee for the critical comments and suggestions 
which have helped in refining the interpretation of our observations and also significantly 
improved the paper. DAR thanks S. Sridhar, R. Nityananda, D. Bhattacharya and K. S. 
Dwarakanath for many useful discussions during the course of the work. We thank C. De 
Pree and D. Balser for useful comments and also for proofreading the manuscript. 
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Table 1. Observational Parameters 


6 cm 3.6 cm 2 cm 0.7 cm 


Date of observations 
Field center RA (J2000) 

Field center DEC (J2000) 

RLs observed - IF1 

- IF2 

Center freq (GHz) - IF1 

- IF2 

Bandwidth (MHz) - IF1 

- IF2 

Velocity range (km s -1 ) - IF1 

- IF2 

Channel separation (km s -1 ) 

- IF1 

- IF2 

Phase calibrator 
Bandpass calibrator 

Flux calibrator 

On-source observing time (hrs) 
Synthesized beam (arcsec) 
Position angle of 
synthesized beam (deg) 

Largest angular size (arernin) 
RMS noise in the spectral 
cube (mJy/beam) 

RMS noise in the 
continuum images (mJy/beam) 


21-OCT-2001 

21-OCT-2001 

19 h 01 m 47 s .l 

19 ft 01 m 47M 

+0TT3'00" 

+01°13'00" 

CllOa 

C92a 

Cilia 

4.8754 

8.3121 

4.7458 

0.78 

1.56 

0.78 

48 

56 

49 

0.75 

0.88 

0.77 

J1851+005 

J1851+005 

J1229+020 

J1229+020 

3C 286 

3C 286 

1.6 

1.8 

14.3 x 13.0 

8.4 x 7.6 

1 

o 

o 

bo 

-17°.2 

5 

3 

1.9 

i.r 

2.0 

1.4 


13-OCT-2001 

13/17-AUG-2004 

19 h 01 m 46 s .4 

19 h 01 m 46 s .4 

+0TT3'24" 

+01°13'24" 

C76a, He76a 

C53a 

14.6927 

42.9638 

3.13 

3.13 

64 

42 


1.0 

1.4 

J1851+005 

J1851+005 

J1924—292 

J1733—130 

J1733—130 


3C 286 

3C 286 

4.0 

2 x 3.3 

4.9 x 4.5 

2.1 x 2.0 

17°.7 

—89°.9 

1.5 

0.7 

1.9 

1.4 

2.9 

0.7 


a Spectral RMS obtained from an off source region. The spectra toward bright continuum emission 
are affected by baseline ripple and hence have higher RMS. 







Tabic 2. 


Parameters of G35.20—1.74 obtained from the continuum images 


Obs. A 
(cm) 

Angular Size 1 
(") 

Flux density 2 

(Jy) 

0.7 

3.7 

2.21 (0.02) 

2 

3.6 

2.49 (0.02) 

4 

< 7.4 

2.41 (0.02) 

6 

< 15.7 

1.95 (0.03) 


1 Deconvolved angular size obtained 
using the AIPS task JMFIT. 

2 The quoted uncertainties are the 
RMS obtained from the residual after 
removing the continuum model. 







Table 3. 


Physical properties of G35.20—1.74 


Distance 

(kpc) 

size 

(pc) 

Te 

(K) 

EM (x 10 7 ) 
(pc cnr 6 ) 

n e (x 10 4 ) 
(cm -3 ) 

U 

(pc cm -2 ) 

Log(iV c ) 

(s- 1 ) 

Spectral Type 1 

3.2 

0.1 

9900 (1400) 

6.4 (0.3) 

2.6 

42.4 

48.4 

08 - 07.5 


1 Spectral type is estimated by assuming a single star is embedded in the UCHii region (Panagia 
1973). 
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Table 4. Parameters of the observed recombination lines 


Obs. A 

Line 

Sl 

AV 

Vlsr 

V f 

v res 

Effective Beanr a 

(cm) 


(mJy) 

(km s -1 ) 

(km s -1 ) 

(km s^ 1 ) 

("x ") 


Toward G35.20-1.72 a 

0.7 

C53a 

14.7 (5.1) 

5.5(2.2) 

42.0(0.9) 

1.4 

6.3 x 6.3 


X53a 9 

80.0 ( 6 . 8 ) 

3.1(0.3) 

43.4(0.1) 

— 

— 

2 

C76 a 

10.7 (1.6) 

4.5(0. 8 ) 

41.9(0.4) 

1.0 

6.3 x 6.3 


He76a 

26.9(0.8) 

21.2(1.3) 

42.0(0.5) 

— 

— 

3.6 

C92a 

(5.2) b 



0.9 

8.4 x 7.6 (—17°.26) 

6 

CllOa & 

( 2 .l) c 



0.8 

14.4 x 13.0 (—0°.33) 


Cilia 






Toward W48A 

2 

He76a 

1.4 (0.2) 

14.8(1.1) 

47.4(0.5) 

1.0 


3.6 


(0.3) 



0.9 


6 


( 1 . 0 ) e 



0.8 



a For 0.7 and 2 cm data, the line parameters are obtained from the spectra averaged 
over the effective beam area centered near the continuum peak. The coordinates of the 
continuum peak are RA: 19 h 01 m 46 s .4, DEC: TORTS* 23”.9, J2000. 

b RMS obtained from the spectrum toward the UCHu region. This spectrum is affected 
by a baseline ripple and therefore the estimated RMS is higher than that obtained from an 
off source position. 

C RMS of the average of Cl 10a and Cilia spectra 

e RMS from the Cilia spectrum 

f Channel separation of the spectra in km s -1 

g Presumably a heavy ion (ion massive than C + ) recombination line; the LSR velocity is 
obtained by assuming an infinite mass for the heavy ion. 
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Table 5. Physical properties of the PDR toward G35.20—1.74 


T PDR 

y.PDR 

1A e 

l n H 

P PDR 

(K) 

(cm” 3 ) 

(x 10” 4 pc) (x 10 6 cm” 3 ) 

(x 10”' dyne cm” 2 



Models with PDR slabs in front and back of the UCHll region 


150 

1500 

- 2000 

0.8 - 

0.5 

5.1 - 

6.8 

5.3 

- 7.1 

200 

1500 

- 2500 

1.2 - 

0.6 

5.1 - 

■ 8.6 

5.6 

- 9.4 

300 

1500 

- 3000 

2.6 - 

0.8 

5.1 - 

10.3 

6.2 

- 12.5 

500 

1500 

- 4000 

5.3 - 

1.0 

5.1 - 

13.7 

7.5 - 

- 20.0 

1000 

1500 

- 4000 

17.0 - 

2.9 

5.1 - 

13.7 

10.6 

- 28.4 


Models with a PDR slab in front of the UCHu region 


200 

2000 

- 2500 

1.3 - 

1.0 

6.8 - 

- 8.6 

7.5 

- 9.4 

300 

2000 

- 4000 

2.2 - 

1.0 

6.8 - 

13.7 

8.3 - 

- 16.6 

500 

2000 

- 5500 

4.5 - 

1.0 

6.8 - 

18.8 

10.0 

- 27.5 

1000 

2000 

- 6000 

14.0 - 

2.6 

6.8 - 

20.5 

14.2 

- 42.6 









Table 6 


LSR velocity of high-density tracers observed toward G35.20—1.74 


Molecular 

transition 

Vlsr 
(km s _1 ) 

f) a 
u res 

(") 

References 

13 CO(l 0) 

45.0 

25 

Churchwell et al. (1992) 

13 CO(2 1) 

46.1 

12 

Churchwell et al. (1992) 

NH 3 (1,1) 

44.5 

40 

Churchwell et al. (1990) 

NH 3 (2,2) 

44.3 

40 

Churchwell et al. (1990) 

CS (2 -> 1 ) 

42.8 

25 

Churchwell et al. (1992) 

CS (5 -f 4) 

43.7 

12 

Churchwell et al. (1992) 


a Angular resolution of the observations 
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19 01 54 52 50 48 46 44 

RIGHT ASCENSION (J2000) 



Fig. 1.— Continuum images of W48A (extended source) and the UCH II region G35.20—1.74 
(compact source) made with the VLA. (a) 6 cm image with an angular resolution of 14".3 
x 13".0 (—0°.8). The contour levels are ( — 1, 2, 5, 10, 20, 32, 64, 128, 256) x 10 mJy/beam 
and the grey scale ranges from —0.01 to 2.55 Jy/beam. The RMS noise in the image is 
2.0 mJy/beam. (b) 3.6 cm image with an angular resolution of 8".4 x 7".6 (—17°.2). The 
contour levels are (—1, 2, 5, 10, 20, 32, 64, 128, 256) x 5 mJy/beam and the grey scale 
ranges from —0.01 to 2.29 Jy/beam. The RMS noise in the image is 1.4 mJy/beam. (c) 2 
cm image with an angular resolution of 4".9 x 4".5 (17°.7). The contour levels are ( — 1, 2, 
5, 10, 20, 32, 64, 128, 256) x 10 mJy/beam and the grey scale ranges from —0.01 to 1.67 
Jy/beam. The RMS noise in the image is 2.9 mJy/beam. (d) 0.7 cm image with an angular 
resolution of 2".l x 2".0 (—89°.9). The contour levels are ( — 1, 2, 5, 10, 20, 32, 64, 128, 256) 
x 3 mJy/beam and the grey scale ranges from —4 to 630 mJy/beam. The RMS noise in the 
image is 0.7 mJy/beam. 




Frequency (GHz) 


Fig. 2.— Flux density vs frequency of the UCHll region G35.20—1.74. The solid curve 
shows the model flux density as a function of frequency for T e = 9900/1 and EM = 6.4 x 10' 
pc cm -6 . The measured flux density at the four frequencies along with ±3<r error bars are 
also marked. 
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Fig. 3.— Spectra of the 53a: and 76a recombination line transitions observed near 2 (bottom) 
and 0.7 (top) cm toward G35.20—1.74. The spectra are obtained by averaging over a 6 ".3 x 
6 ".3 region near the continuum peak (RA = 19 ft 0l m 46 s .4, DEC = +01 o 13 , 24 ,, , J2000) in the 
2 and 0.7 cm images. The line features near +42 km s -1 are the carbon lines. The He76a 
with an LSR velocity of +42.0+0.5 km s -1 is the second feature in the spectrum shown on 
the bottom panel. Only part of the helium line profile is detected due to limited bandwidth 
of the observation. The line feature near +30 km s^ 1 in the top spectrum is presumably a 
heavy ion (ion massive that C + ) recombination line. The LSR velocity is with respect to 
the carbon lines. The result of a two component Gaussian fit to the line features is shown 
by the solid line. 
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Fig. 4.— Carbon line flux density as a function of frequency for class (a) and (b) models (see 
§4.2) are shown on Figs, (a) and (b) respectively. The solid and dashed curve corresponds 
to gas temperatures of 500 and 200 K respectively. The electron density used to obtain all 
curves are 2500 cm -3 . The upper limits at 6 and 3.6 cm wavelengths and the line flux 
densities of C53a and C76 a transitions with ± la error are also marked. At frequencies 
below ~ 10 GHz, the 200 K model predicts the carbon RLs in absorption. The parameters 
of such models are constrained by using the 6 cm spectral uncertainty (3cr) as a lower limit 
for the model flux densities at this wavelength, which is also shown in the figure. At 3.6 cm 
the upper limit obtained from off source spectra (3.3 mJy) is shown with a circle. Note that 
the 200 K model flux densities are consistent with this upper limit. 
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Observer 


Fig. 5.— Schematic of the UCHll region G35.20—1.74, the associated PDR and the molec¬ 
ular cloud. The hydrogen RL at 14.6 GHz originates from the UCHll region. The central 
velocity of the H76cr transition (+47.9±1.2 km s _1 ) indicates the mean velocity of the ion¬ 
ized gas with respect to the LSR. The carbon RL observed at 14.6 GHz is detected from 
the region of the PDR marked by the curving line, since the line emission is dominated by 
stimulated emission. The central velocity of C76a line (+41.9±0.4 km s _1 ). The mean LSR 
velocity of the observed molecular lines is +44.4±1.0 km s -1 . Thus relative to the molecular 
cloud the PDR is moving into the cloud at —2.5±1.1 km s -1 and the ionized gas is moving 
at +3.5±1.6 km s -1 . 



